home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / COMQR.z / COMQR
Text File  |  1996-03-14  |  4KB  |  67 lines
  1.  
  2.  
  3.  
  4. ____CCCCOOOOMMMMQQQQRRRR((((3333FFFF))))                                                          ____CCCCOOOOMMMMQQQQRRRR((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      COMQR, SCOMQR   -  EISPACK routine.  This subroutine finds the
  10.      eigenvalues of a COMPLEX upper Hessenberg matrix by the QR method.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  ccccoooommmmqqqqrrrr((((nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, hhhhrrrr,,,, hhhhiiii,,,, wwwwrrrr,,,, wwwwiiii,,,, iiiieeeerrrrrrrr))))
  15.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, iiiieeeerrrrrrrr
  16.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn hhhhrrrr((((nnnnmmmm,,,,nnnn)))),,,, hhhhiiii((((nnnnmmmm,,,,nnnn)))),,,, wwwwrrrr((((nnnn)))),,,, wwwwiiii((((nnnn))))
  17.  
  18.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssccccoooommmmqqqqrrrr((((nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, hhhhrrrr,,,, hhhhiiii,,,, wwwwrrrr,,,, wwwwiiii,,,, iiiieeeerrrrrrrr))))
  19.           iiiinnnntttteeeeggggeeeerrrr          nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, iiiieeeerrrrrrrr
  20.           rrrreeeeaaaallll             hhhhrrrr((((nnnnmmmm,,,,nnnn)))),,,, hhhhiiii((((nnnnmmmm,,,,nnnn)))),,,, wwwwrrrr((((nnnn)))),,,, wwwwiiii((((nnnn))))
  21.  
  22.  
  23.  
  24. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  25.      On INPUT
  26.      NNNNMMMM must be set to the row dimension of two-dimensional array parameters
  27.      as declared in the calling program dimension statement.
  28.  
  29.      NNNN is the order of the matrix.
  30.  
  31.      LLLLOOOOWWWW and IGH are integers determined by the balancing subroutine  CBAL.
  32.      If  CBAL  has not been used, set LOW=1, IGH=N.
  33.  
  34.      HHHHRRRR and HI contain the real and imaginary parts, respectively, of the
  35.      complex upper Hessenberg matrix.  Their lower triangles below the
  36.      subdiagonal contain information about the unitary transformations used in
  37.      the reduction by  CORTH, if performed.  On OUTPUT The upper Hessenberg
  38.      portions of HR and HI have been destroyed.  Therefore, they must be saved
  39.      before calling  COMQR  if subsequent calculation of eigenvectors is to be
  40.      performed.
  41.  
  42.      WWWWRRRR and WI contain the real and imaginary parts, respectively, of the
  43.      eigenvalues.  If an error exit is made, the eigenvalues should be correct
  44.      for indices IERR+1,...,N.
  45.  
  46.      IIIIEEEERRRRRRRR is set to ZERO       for normal return, J          if the J-th
  47.      eigenvalue has not been
  48.         determined after a total of 30*N iterations.  Calls CSROOT for complex
  49.      square root.  Calls PYTHAG(A,B) for sqrt(A**2 + B**2).  Calls CDIV for
  50.      complex division.  Questions and comments should be directed to B. S.
  51.      Garbow, APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.